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Abstract. In this article we describe our experiences with a parallel SIN- 
GULAR - implementation of the signature of a surface singularity defined by 
2^ + 3(x,y) = 0. 

1. Introduction 

Let g G C[a;, y] define an isolated curve singularity at G and / := +g{x, y) 
for N > 2. The zero-set V :— V{f) C C'^ of / has an isolated singularity at 0. 
For a small e > let K := V{f - e) C be the Milnor fibre of {V, 0) and 
s : H2(ye,R) X ff2(K,K) — ^ K be the intersection form (cf. [1], [10], [11], [18]). 
H2{Vi;,R) is a /z-dimensional R-vector space, /x the Milnor number of {V,0) (cf. 
[Tj, [7], [H], [5]), and s is a symmetric bilinear form. Let a{f) be the signature of 
s, called the signature of the surface singularity {V,0). Formulas to compute the 
signature are given by Nemethi (cf. [T3], [Hj) and van Doorn, Steenbrink (cf. [1]). 
We have implemented three approaches in Singular (cf. [3], [6]) using Puiseux 
expansions, the resolution of singularities and the spectral pairs of the singularity. 
In section [2] we will explain the different formulas to compute the signature, and 
finally we give examples and timings of our implementation in section |3l 

2. The signature of {V, 0) in terms of N and g 

The proofs of the following formulas (cf . Propositions 12. 1[ 12. 2[ 12.41 and 12. 7p can 
be found in the corresponding papers of Nemethi (cf. [T3], |14)). 

2.1. Approach 1: Puiseux Pairs. For the first approach assume that {V{g), 0) C 
(C^, 0) is the germ of an irreducible curve singularity. Let (mi, ni), . . . , (m^, ni) be 
the Puiseux pairs of g and define a sequence {aj}i=i_...,£ by 

ai = mi and a,+i = m^+i - n^+i • (m; - n, ■ a,). 

Moreover, we set dg — 1 and di — gcd(7V, n^+i • • • rtg) ior 1 < i < i. 

Proposition 2.1. a{f) = Li ^« ' {^"'^ + y"* + z^/'^') . 

The signature of a Brieskorn polynomial x''^ + y'^^ + z'^^ can be computed com- 
binatorially. Let St = #{(A:i, fcz, ^a) S | 1 < < aj ~ l,t < X^^^i < ^ + 1} 
for t£No. 

Proposition 2.2. cr(x'=i + + z''') ^ So - Si + 82- 
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Remark 2.3. The Singular implementation of the first approach bases on the 
procedure invariants (cf. library hnoether . lib) to obtain the Puiseux pairs of 
g. The command list L = invariants (g) provides a list L in which the third 
and fourth entry contain the necessary data for our application. The required 
combinatorics of Proposition 12. II resp. 12.21 had to be implemented in spring 2011. 

2.2. Approach 2: Resolution. In the following 0) doesn't need to be irre- 
ducible. For the second approach we use the resolution of the singularity {V{g), 0). 
Let V = W ]J ^ be the vertices of the resolution graph, W the vertices correspond- 
ing to the exceptional divisors and A the vertices corresponding to the resolved 
branches. Let £ = {{vi,V2) \ vi,V2 G V} be the set of edges of the resolution 
graph. Let {TO^jugyy be the sequence of total multiplicities and set rua = 1 if 
a G A. For it; e W let = gcd(mi, \ v G Vw U {w}) and for e = (wi, U2) £ £ let 
uie = gcd(mi,i, m^jj. 

Proposition 2.4. a{z^ + g) = 77(5, N) - N ■ r]{g, 1) and^ 

ees wew 

+4- E E E((^))-((^))- 

w has more {v,w)^£ 
than 2 edges 

Moreover, it holds Y.ee£ (scdiK^me) - l) = Y.,,^^ {gcd{K, M^) - l) tf{V{g),0) 
is irreducible. 

Remark 2.5. The Singular implementation of the second approach bases on 
the procedure totalmultiplicities (cf. library alexpoly . lib) to obtain nec- 
essary information about the resolution of {V{g),0). The command list L = 
totalmultiplicities (g) ; provides a list L of the following format: L[l]: in- 
cidence matrix of the resolution graph, L [2] : list of the sequences of the total mul- 
tiplicities corresponding to the branches, L [3] : list of the multiplicity sequences of 
the branches. The required combinatorics of Proposition l^^ had to be implemented 
in spring 2011. 

2.3. Approach 3: Spectral Pairs. The third approach uses the spectral pair^ 
of the singularity {V{g),Q). Therefore let 

Spp{g) = E ■ 

{oi,w) 

represent the spectral pairs where Sq = if a ^ Z and Sq = 1 if a G Z, and 
Sp{9) = E(a,^06Spp(g)(") be the spectrum of g. 



^The command list L = invariants (g) ; returns a list L of the following format: L[l]: char- 
acteristic exponents, L[2]: generators of the semigroup, L[3]: first components of Puiseux pairs, 
L[4]: second components of Puiseux pairs, L[5]: degree of the conductor, L[6]: sequence of 
multiplicities. 

•^For a; G Q we denote by {x} the fractional part of x and 

{x}-\, ifx^Z 

0, if X e z ■ 



((^■)) 



The definition of the eta-invariant ri{f, K) can be found in I12| . 
■^A definition of spectral pairs can be found in [9]. 



Remark 2.6. Note that a is a spectral number, i.e. exp(— 27ria) is an eigenvalue 
of the monodromy. {V{g),0) is reducible if and only if is a spectral number (cf. 
[1]). The spectral numbers are situated in the interval (—1, 1) and the spectrum is 
symmetric (a is a spectral number if and only if —a is a spectral number, cf. |,9,). 
If the Newton polygon of g is non-degenerate the spectral pairs can be computed 
combinatorially using the Newton polygon (cf. [16], [E]). There is a formula to 
compute the spectral pairs from the data of the resolution (cf. [17]). 

Proposition 2.7. a{z^ + g) — v{9i ^) ~ ^ " v{9j 1) '^'^'^ 



a=iO,KaeZ a>0,Ka^Z 
(a,2)£Spp(g) (a,w)£Spp(g) 

Remark 2.8. The Singular implementation of the third approach bases on the pro- 
cedure sppairs (cf. library gmssing . lib) to obtain the spectral pairs of {V{g),Q). 
The command list L = sppairs (g) ; provides a list L of the following format: 
L[l]: set of spectral numbers {ai,...,ar}, L[2]: set of weights {wi, . . . , w^}, 
L [3] : set of multiplicities {hi, . . . , h^} such that Spp{g) = ' i'^j^ ^j)- The 

required combinatorics of Proposition 12 . 71 had to be implemented in spring 2011. 

2.4. Theoretical Comparison. The topological type of a plane curve singularity 
defined by g{x, y) = can be described by the Puiseux pairs (Approach 1) of the 
branches and their intersection multiplicities or equivalently by discrete invariants 
of the resolution (Approach 2) There are combinatorial formulas to get from one 
description to another (cf. [7]). Moreover, the spectral pairs (Approach 3) which are 
topological invariants introduced by Arnold (cf. [T]) and Steenbrink (cf. [T^) can 
also be computed combinatorially from the resolution data (cf. [I?]). Consequently, 
all of the three approaches to compute the signature of the surface + g{x, y) = 
as described above are based on the knowledge of three different finite sets of invari- 
ants which are related in a combinatorial way. The essential difference concerning 
these approaches is the method to compute the set of the corresponding invariants. 

The spectral pairs (Approach 3) can be computed from the mixed Hodge struc- 
ture. This requires several standard basis computations of certain modules over lo- 
cal rings which is the bottleneck of this approach. The Singular library 
gmssing. lib is designed for computing the mixed Hodge structure for hypersurface 
singularities of any dimension. This is one reason why computing the spectral pairs 
using this method is usually comparatively slow. 

The Puiseux pairs (Approach 1) of the branches and their intersection multiplic- 
ities resp. the resolution graph (Approach 2) and the multiplicity sequence can be 
computed via Hamburger-Noether expansion (cf. [1]) or resolution of the curve sin- 
gularities. Hence, both approaches are similarly time-consuming. They only need 
Grobner basis computations if field extensions of Q are necessary to compute the 
Hamburger-Noether expansion resp. the resolution. Anyway, the field extensions 
are the bottlenecks of these approaches. 
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3. Examples and timings 



In this section we provide examples on which we time the three approaches 
as described in section [5] to compute the signature of a surface singularity + 
g{x,y) — 0. The corresponding procedures are implemented in Singular in the 
library surf acesignature . lib. Timings are conducted by using Singular 3-1-3 
on an Intel® Xeon® X5460 with 4 CPUs, 3.16 GHz each, 64 GB RAM under the 
Gentoo Linux operating system. 

Example 3.1. We consider the following polynomials: 

gi = x^^ ~ 21x^^ + Sx^^y - 6a:" - IGx^^y + 20x^^y'^ - x^'^ + 8a:"?; - 36x^°y'^ 
+ 24:x^y^ + 4x^2 - 16a:^y^ -I- 26a;^y'' - 6a;^?;'* + 8x^y^ + ix^y^ - y^, 



92 


= 9U 


x^'y'^ 


93 


- (y' 


+ 2x^y 


94 




a:20y20 


95 






96 


= x^"^ 





The curve singularities in (C^, 0) defined by gi resp. g2 resp. 173 are analytically 
irreducible with Puiseux pairs (3, 2), (7, 2), (15, 2) resp. (3, 2), (7, 2), (15, 2), (67, 3) 
resp. (3, 2), (7, 2), (113, 3). 

The curve singularities in (C^, 0) defined by 174 resp. 175 resp. ge are analytically 
reducible since they are intersections of 40 resp. 10 resp. 20 lines at the origin. 
Consequently, the first approach is not applicable for these examples. Furthermore, 
the polynomials 54,35,55 are defined over Q, whereas the resolution is only defined 
in field extensions of degree 8 resp. 10 resp. 20 over Q. 

Computations reveal the following results and corresponding timings which are 
summarized in Table [1] The symbol " > 14 h" indicates that the computation did 
not terminate after more than 14 hours. All timings are given in seconds (s). 



N 


9i 


(^{z^ +9r{x,y)) 


Approach 1 [s] 


Approach 2 [s] 


Approach 3 [s] 


5 


91 


-168 








> 14 h 


5 


92 


-1620 


174 


183 


> 14h 


6 


93 


-940 


2908 


2912 


> 14h 


5 


94 


-3192 




19 


> 14h 


6 


95 


-189 




22 





6 


9e 


-779 




14542 
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Table 1. Results and total running times for computing the sig- 
nature of the surface singularity given by the considered examples 
(cf. Example 13. 1|) via all approaches as described in section [2] 



In addition, we summarize the maximal memory allocated from system during 
the considered computations in Table [2] The memory consumption is given in 
Megabyte (MB). 
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N 


9i 


Approach 1 [MB] 


Approach 2 [MB] 


Approach 3 [MB] 


5 


91 


1 


1 


> 300 


5 


92 


3442 


3442 


> 1100 


6 


93 


7723 


7723 


> 1100 


5 


54 




33 


> 2300 


6 


55 




8 


4 


6 


56 




156 


72 



Table 2. Maximal memory allocated from system while comput- 
ing the signature of the surface singularity given by the considered 
examples (cf. Example 13. 1|) via all approaches as described in sec- 
tion El 



Remark 3.2 (Algorithmic Conclusion). Our experiments reveal that there exist 
examples where Approach 1 and Approach 2 are almost equivalent regarding time 
consumption resp. memory allocation, but Approach 3 is much slower (cf. gi, g2, 
g^). Furthermore, there exist examples where Approach 1 is not applicable, but 
Approach 2 consumes more time resp. allocates more memory than Approach 3 
(cf. (74), and vice versa (cf. 175, gg). Consequently, it is reasonable to summarize 
all approaches in one algorithm which computes the signature via every approach, 
if possible, in parallel such that the fastest approach wins and returns the result. 
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